Enhancement of the upper critical field by nonmagnetic impurities in dirty two-gap 

superconductors 



(N 
O 

o 
o 

Q 
in 



o 
o 

u 



O 
o 



> 

O 



I 

o 
o 



A. Gurevich 

Applied Superconductivity Center, University of Wisconsin, Madison, 

(Dated: February 1, 2008) 



Wisconsin 53706 



Quasiclassic Uzadel equations for two-band superconductors in the dirty limit with the account of 
both intraband and interband scattering by nonmagnetic impurities are derived for any anisotropic 
Fermi surface. From these equations the Ginzburg-Landau equations, and the critical temperature 
Tc are obtained. An equation for the upper critical field, which determines both the temperature 
dependence of Hc2{T) and the orientational dependence of Hc2(d) as a function of the angle 9 
between H and the c-ajcis is obtained. It is shown that the shape of the Hc2{T) curve essentially 
depends on the ratio of the intraband electron diffusivities D\ and D\, and can be very different 
from the standard one-gap dirty limit theory. In particular, the value HciiQ) can considerably 
exceed 0.7TcdHc2/dTc, which can have important consequences for applications of MgB2- A scaling 
relation is proposed which enables one to obtain the angular dependence of Hc2 (Q) from the equation 
for Hc2 at H||c. It is shown that, depending on the relation between D\ and D2, the ratio of the 
upper critical field H^},2/H^ for H||a6 and H _L a6 can both increase and decrease as the temperature 
decreases. Implications of the obtained results for MgB2 are discussed. 

PACS numbers: PACS numbers: 74.20.De, 74. 20. Hi, 74.60.-w 



I. INTRODUCTION 



Recent discovery of MgB2 has renewed interest in su- 
perconductors with multicomponent order parameters. 
In particular, extensive ab-initio calculationsaB, along 
with ample experimental jfividences based on STMBq, 
point contactO, and Ramantl spectroscopy, heat capacityCl 
and rf responsetl measurements strongly indicate that 
MgB2 has two distinct s-wave superconducting gaps 
A<,(0) w 7.2mV and A^(0) w 2.3mV. These gaps re- 
side on different disconnected sheets of the Fermi sur- 
face (FS), which comprises nearly cylindrical 2D parts 
formed by in-plane cr antibonding p^y orbitals of B, and 
a more isotropic 3D tubular network formed by out-of- 
plane tt bonding and antibonding orbitals of B. So 
far all attempts to increase the critical temperature Tc of 
MgB2 by doping have been unsuccessfula. On the othet. 
hand, the significant potential of MgB^ for applicationsE2l 
is limited by comparatively low upper critical fields 
ff^-L(O) ~ 3 - 5T and pS'Uo^p^i^p^mXipf.iie^ 
MgB2 single crystalfJliOliallil'BByyyMEil, where 
the indices _L and || correspond to the magnetic field H 
perpendicular and parallel to the ab plane, respectively. 
As far as Hc2 is concernedj-i|t can be increased by adding 
nonmagnetic impuritiesc3'E3, which are especially effec- 
tive in the dirty limit, 2'KkBTc < ti/r, where r is the 
elastic scattering time, h is the Planck constant, and kB 
is the Boltzmann constant. In this case, there is a sim- 
ple universal relation between the zero-temperature value 
-f^c2(0), the slope i?^2 — dlL^/dT at Tc, and the normal 
state residual resistivity p„E3: 



where Np is the density of states at the FS and — e is the 
electron charge. The fact that adding nonmagnetic im- 
purities increases Hc2 has been used to significantly im- 
prove Hc2 of NbTi and NbsSn compoundf£3. The same 
approach has been also taken to improve the high-field 
performance of MgB2 in which the additiorLai scatter- 
ing was introduced by protoncj and ncutroncJ irradia- 
tion, and atomic substitutions on both-JS-and Mg sites in 
thin filmsBS, bulk samplesO and wirescIlEs. Fot-instance, 
in highly resistive c-axis oriented Mgi?2 filmsc3, p„ was 
increased from the nominal value ^ l/ii7cmEj to more 
than 700/^i7cm, which in turn increased the slopes _ffc2 
to rather high values iJ^2_L ^ IT/K and if^2|| ~ 1-8T/K, 
while reducing Tc down to ~ 30K. Based on these num- 
bers, the extrapolation (]l|) would give the upper critical 
field, H^^(D) w 20T, which is stiU below Fc2(0) « 30T 
of iVfoaS'nci However, Eq. (P is a result of the one- 
gap theoryE3 which does not take into account the mul- 
tiple scattering channels in two-gap superconductors. In 
this paper I show that -ffc2(0) in two gap superconduc- 
tors can be significantly higher than what follows from 
Eq. (|l|), thus the upper critical fields of MgB2 alloys 
can in fact exceed Hc2{0) of Nb^Sn even for thp-,yaiues 



Hc2± 



IT/K which have already been achievedo' 



Hc2{0) = 0mTcH'c2 



H'c2 



Aeki 



■NpPn, (1) 



The theory of the twoiband superconductivity has 
been developed long agotil, and more ijec |i j'n|tlj fj^eyral- 
ized to include the effects of impiU|d||ief£3Ba and 
strong electron-phonon couplingpEjO'Lj. The Fermi sur- 
face of MgB2 provides three different impurity scattering 
channels: intraband scattering within each a and tt FS 
sheets and interband scattering between them. Strong 
impurity scattering in the dirty limit causes intraband 
electron diffusion over the respective FS sheets, which re- 
duces the intrinsic anisotropy of Ag- and A^r , but does not 



2 



affect Tc, in accordance with the Anderson theorem. By 
coptpst, the pairbreaking interband scattering reduces 
Tct3H, but its effect in seems to be not verii-pro- 

nounccd due to orthogonaHty of a and tt orbitalso. It 
is the multiple scattering channels, which provide the es- 
sential flexibility to increase the Hc2 of A/.gi?2 to a much 
greater extent than in one-gap superconductors not only 
by the usual increase of p„, but also by optimizing the 
relative weight of the a and tt scattering rates. The latter 
could be done by selective atomic substitution on B sites 
(with C, O, etc. ) which mostly affect the inriplane a 
scattering, and Mg site (Mg vacancies ,pA.l, etc)El, which 
mostly affects out-of-plane tt scatterin^J. 

In this work I use a weak coupling multiband BCS 
modefij to develop a theory of magnetic properties of 
two- gap superconductors in the dirty limit, taking MgB2 
as an example. This approach generalizes the well-knoBm 
theory developed for dirty one-gap superconductorscj. 
The paper is organized as follows. In section II the qua- 
siclassic Uzadel equations are obtained for an anisotropic 
two-gap superconductor with the account of both inter- 
band and intraband scattering by nonmagnetic impuri- 
ties. In section III these equations are used to derive the 
Ginzburg-Landau equations and the critical temperature 
Tc in the dirty limit. In section IV the linearized Uzadel 
equations are solved exactly to obtain a general equation 
for the upper critical field parallel to the c-axis. It is 
shown that the temperature dependence n£ Hc2{T) can 
be very different from the one-band theorycj, and Hc2 (0) 
can be considerably higher than what follows from Eq. 
(||). In section V I calculate the angular dependence of 
Hc2{d) and propose a scaling rule, which shows that the 

anisotropy ratio H^c2l -^02 can both increase and decrease 
with temperature, depending on the relation between in- 
traband diffusivities. In section VI the results of this 
work are used to address the observed temperature and 
orientational dependences of Hc2 of MgB2. 



II. UZADEL EQUATIONS 

A powerful tool for investigating inhomogeneous 
states of supercatuiuctors is the quasiclassic Eilen- 
berger equationsEjCj for the Green functions /(k, r,^), 
/+(k, r,^), ^(k, r,ijj) and the order parameters A(k, r) 
which depend on the coordinates r, the Matsubara fre- 
quency Lo = 7rT(2n -I- 1), and the wave vector k on the 
FS (hereafter the units with = ti = I are used). The 
essential dependence of /(k, r,(jj), ^(k, r, oj) and A(k, r) 
on the direction of k, makes them sensitive to the shape 
of the FS, which greatly complicates solving the nonlin- 
ear Eilenberger equations. Thus, the Eilenberger equa- 
tions are not so useful in the case of complicated FS, 
such as the tubular network of disconnected FS sheets of 
MgB2 for which the numerical solutions of the Eliash- 
berg equation have shown the significance of the FS 
anisotropycl. However the situation simplifies in the dirty 
limit for which the impurity scattering within each FS 



sheet averages out the angular dependences of /(k, r,ti;), 
/^(k, r, u), g(\i, r, w), making them independent of k. In 
this case the Eilenberger equations reduce to much sim- 
pler Uzadel equation^ in which all microscopic details 
are hidden in the electronic diffusivity tensors D"f for 
each m-th FS sheet and the interband scattering rates 
7mm' which reflect the underlying FS symmetryuS. 

In this paper I consider a dirty two-gap anisotropic 
superconductor in the simplest case of two disconnected 
sheets 1 and 2 of the FS._This model corresponds to the a 
and TT bands in Mg_B2cl0, for which the superconducting 
gaps take the constant values Ai and A2 on the sheets 
1 and 2, respectively. The Uzadel equations derived in 
Appendix A have the form 



2Aigi + 7i2(.9i/2 - 32/1) 
2Luf2 - Df[g2llcJVpf2 - f2Vo.Vpg2] 
2A2.g2 + 72i(.92/i - .91/2), 



(2) 



(3) 



Eqs. (H) and (||) are supplemented by the self-consistency 
equations for the order parameters A^ = |Am| exp(i(y9m): 



Am — 2ttT ^ ^ ^ ^ ^vam' fm' (l". 



UJ 



LJ>0 m 

normalization condition 



\fr, 



9m 



and the expression for the supercurrent density 



(4) 



(5) 



(6) 



Here the band index m runs from 1 and 2, the functions 
/„j(r, uj) and gmi^, in the m-th band depend on r and 
UJ but not on k, D"J^ are the intraband diffusivity tensors 
due to nonmagnetic impurity scattering, 7mm' are the 
interband scattering rates, iVm is the partial density of 
states, n = V -I- 27riA/0o, A is the vector potential, (j)o 
is the flux quantum, and the summation over the Greek 
Cartesian indices is implied. 

Eqs. (^) contains the matrix of the BCS supercon- 
ducting coupling constants Amm' = ^mm' ~ A^mm' , where 



•^mm' electron-phonon constants, and /imm' is the ma- 
trix of the Coulomb pseudopotential. Here the diagonal 
terms An and A22 quantify the intraband superconduct- 
ing coupling, and off-diagonal terms A12 and A21 describe 
the interband coupling. The eigenvalues of Amm' are as- 
sumed positive, the band 1 having the highest coupling 
constant An. The indices 1 and 2 thus correspond to 
the a and tt bands of MgB2 for which ab-initio calcula- 
tions yield X^^ w 0.81, X^ttt ~ 0.285, Ao-tt ~ 0.119, and 
Ajro- ~ 0.09123. The mixedj-cpmponents A12 and A21 sat- 
isfy the symmetry relatiorH (see also Eq. (A23)): 



TViAi 



(7) 
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where A^i and N2 are partial densities of states in the 
bands 1 and 2 (iV^ « l.iN„ for MgB2)- A similar ap- 
proach based on the Uzadel equations was recently ppa- 
posed to describe vortices in two-gap superconductorsc3. 

Formulas for the interband scattering rates 7mm' and 
the intraband diffusivity tensors expressed in terms 
of microscopic material parameters are given in Appendix 
A. They will be used in the next sections as input param- 
eters, which can be either calculated from first principles 
or extracted from the observed temperature dependences 
of Hc2{T) and resistivity p, as discussed below. Here we 
just emphasize the fact that, although both D'^^ and 
reflect the underlying symmetry of the FS, the fea- 
tures of atomic orbitals which form the bands 1 and 2 
can manifest themselves in very different properties of 
D"^ and £'2'^. For instance, the principal value D^f' 
along the c-axis in MgB2 is much smaller than two prin- 
cipal values ZJ^"' and D^J''' in the ab plane, because of the 
nearly 2D nature of the cr-band formed by in-plane bond- 
ing asd antibonding p^y orbitals of B, as was shown by 
STMB. By contrast, the difference in the principal values 



and dI^^ is less pronounced for the 3D 7r-band 



formed by the out-of-plane Pz orbitals of B. The result- 



ing relation Da /Da ^ Di^' /Dl^' can manifest itself 
in anomalous temperature dependence of the anisotropy 
of i?c2, as shown in the next sections. 



term in the right hand side of Eq. 
summing up over uj: 



I) into Eq. (|) and 



(12) 



where the asterisk means complex conjugation. As fol- 
lows from Eq. (|l^), the current density near Tc is a sum 
of independent intraband contributions. 

Equilibrium equations for the order parame- 

ters Am are not yet the GL equations which are obtained 
by varying the free energy functional, 5F/dA* = 0. The 
reason is that the off-diagonal terms with respect to the 
band index m in Eqs. (|)-® violate the necessary 
symmetry conditions which the variational equations 
5F/5A* — must satisfy. To find the proper linear com- 
binations of Eqs. and (|l^) which satisfy all necessary 
symmetry conditions, we first solve Eqs. for the 

quantities A^i-Ri = (A22A1 — Xi2A2)Ni/w and A^2^2 = 
(All A2 — A2iAi)iV2/w, where w — A11A22 — A12A21 is the 
determinant of the matrix Amm'- Expressing then Rm 
via Am using Eq. (pi]), we obtain the GL equations for 
an anisotropic two-band superconductor: 

- aiAi + 6iAi|Ai|2 + cf n„n^Ai + 5A2 = 0, (13) 
-a2A2 -1- 62A2IA2I2 + cf n„n^A2 + 5A2 = 0. (14) 

Here the GL expansion coefficients are given by 



III. CRITICAL TEMPERATURE AND 
GINZBURG-LANDAU EQUATIONS 

The necessity to satisfy the self-consistency condition 
(^ complicates solving the nonlinear Eqs. (||)-(||) for in- 
homogeneous Am(r)- The situation simplifies near Tc 
where Eqs. (|^)-(P reduce to the Ginzburg-Landau (GL) 
equations for the order parameters Am- We derive the 
GL equations neglecting the interband scattering terms 
712 T^Tc which are usually small \fL.MgB2 due to or- 
thogonality of the a and vr orbitalsEJ. Expanding the 
solution of Eqs. (^) and (||) in powers of Am and its spa- 
tial derivatives, and using — ^ — |/mP/2, we obtain 



ai 



0.2 



In 



In 



ttT 

2-/LUD 



a/3 



ttT 



22 
W 

An 

w 



Ni, 61 



No 



ttNiD^''/8T, 



af3 _ 



7C(3)A^i 
87r2r2 ' 

7C(3)A^2 
87,27-2 ' 

ttN2D^^/8T, 



b2 = 



g = X21N2/W = X12N1/W. 



(15) 

(16) 

(17) 
(18) 



The equality in Eq. ( p^ ) results from Eq. (^ for the in- 
terband coupling constants A12 = ViN2, and A21 — ViNi. 
The GL equations (13) and dll) can now be obtained by 
varying the free energy F = J [-Fi + F2 + Fi]dV, which 
contains the usual GL intraband contributions 



A D"'^ A A 2 

^m -^m TT TT A ^m|^m| 
Jm h — ^liail/sAm 7^ (o) 



Inserting Eqs. (^) into Eq. (Q) and summing up over ui, 
results in the following equations 



i^m = -^lAmP + ^IAml^ + ^H^AmH^A;,, 



and the interband interaction term 

i^, = ^(AiA; + A2A* 



(19) 
(20) 



Rm — Ajjil 



Ai = Aiii?i -t- Ai2i?2, 

A2 = A2ii?i + X22R2, 

" -^m ^ 

8T 



-n^H/jA 



7C(3) 
87r2r2 



(9) 
(10) 

AmpAm, (11) 



where i?m = J2uj hi, I ^ \n{2^ujD /ttT), In 7 = -0.577 is 
the Euler constant, and C(3) = 1.202. The supercurrent 
density at T — Tc <C Tc is obtained by inserting the first 



Notice that Eq. (|lj) for the current density obtained 
from the Uzadel equations coincides with the general 
expression J = SF/SA, where F is defined by Eqs. 
{0 and (|0|). Static GL equaiions (^-(|ol) and their 
time-dependent generalizatiorO have been suggested 
phenomenologically to describe various distributions of 
the order paranaeter and vortex properties in two-gap 
superconductors! 
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Eqs. (|T^) and (^4|) give the following equation 0102 — 
= for the critical temperature T^, which results in 
the quadratic equation for Ic — hi{2ju!D/T^Tc): 



1 - X+l^ + wll = 



(21) 



where A + = Aii±A22, and w — A11A22 — A12A21. Solution 
of Eq. reproduces the well-known result by Suhl, 

Matthias and WalkeiEl: 



r,o = 1.14cji5exp[-(A+ - Ao)/2w] 



(22) 



where Aq = (A^^ +4Ai2A2i)i/2. As follows from Eq. (pi 



the interband coupling always increases Tc as compared 
to noninteracting bands (A12 — A21 = 0), while the in- 
traband impurity scattering does pot affect T^, in accor- 
dance with the Anderson theorerrO. 



IV. UPPER CRITICAL FIELD FOR H||c 

Now we turn to the calculation of the upper critical 
field Hc2 applied along the c-axis of a hexagonal crystal. 
Here Hr2 is the maximum eigenvalue of the linearized 
Eqs. (||) and (^) for / ^ 1 and g ^ 1. In Appendix 
B, these equations are solved exactly and an equation for 
Hc2 is obtained for arbitrary relation between Di , D2 and 
7mm' ■ This general equation for Hc2 is very cumbersome, 
so to make the essential effects of two-gap superconduc- 
tivity more transparent, I consider here a much simpler 
case of negligible 7m„i' for which the nonmagnetic impu- 
rity scattering does not affect Tc- If 712 = 721 = 0, the 
linearized Eqs. (||) and (||) take the form 



2^/1- 
2c^/2 



2Ai, 
= 2A,. 



(23) 
(24) 



For H||c, the relevant solutions of Eqs. (|2^) and (|24| ) 
depends only on isotropic in-plane diffusivities, 1?^'^ = 
DmSaf3- In the gauge Ay — Hx the solutions of Eqs. (|^), 
(||), and (H) are 

fm{x,Uj) = Arn{x)/{uj+nHDrn/<l)o), (25) 

A^{x) = A,„cxp(-7riJa;V(/.o). (26) 

Inserting Eqs. ( |25| ) and (|2^) into the gap equation (^) 
yields two linear equations for Ai and A2 in which the 
summation over uj is performed using the identity 



27rTy — 



w>0 



uj + n ttT ^(27rr) 



C/(x) = V(l/2 + a;)+^/>(l/2), 



(27) 
(28) 



where ip{x) is the di-gamma function. The equations for 
Ai and A2 become 

Ai = Aii[/ - U{h)]Ai + Xi2[l ~ U{T]h)]A2, (29) 
A2 = A22[^ - 1/(77/1)] A2 + X2i[l - U{h)]Ai, (30) 



where I = ln(27W£,/7rT), h = Hc2Di/2(j)oT\ and 77 = 
D2/D1. The solvability condition of Eqs. ( |29| ) and ( ]30| ) 
gives an equation for Hc2, in which it is convenient to 
express ojd via Tc using Eqs. ( pT]) and (p^. As a result, 
the equation for Hc2 takes the form 



2ao [In t - 
a2[\nt + U{7]h)] 



U{h)][\nt + U{rih)] ■ 
-ai[liit + U{h)] =0. 



(31) 



where ai = 1 + A_/Ao, 02 = 1 — A_/Ao, and oq = w/Xq. 
For equal diffusivities, = 1, Eq. (^ reduces to 
the equation Int^ U{h) = for Hc2 in one-gap dirty 
superconductorsE3. A somewhat different form of Ea. 
( ^ ) was recently obtained from the Uzadel equations^. 

Shown in Fig. 1 are Hc2{T) curves calculated from 
Eq. ( pT| ) for different diffusivity ratios 77 = D2/D1. One 
can clearly see the evolution of the temperature depen- 
dence of Hc2 (T) from the dirty- limit one-band behavior 
at Di — D2, to rather different Hc2{T) curves which have 
portions with both upward and downward curvatures for 
either Di <C D2 and Z?i ^ Z?2- In the latter two cases, 
the zero-temperature value i?co(0) can be significantly 
higher than the one-gap extrapolation (0). To show this, 
we obtain expressions for Hc2{T) near Tc and T = 0. For 
T K, Tc, Eq. ( |3l| ) can be expanded in small terms ~ 
using U{h) « 7r^/2, and then solving for Hc2'- 



Hc2 — 



){Tc-T) 



7r2(ai£)i 021)2)^ 



(32) 



For Di = D2, Eq. (|3g) yields the result of the one- 
band theory, Hc2 = 4(Tc - T)(I)o/tt^ cM. Furthermore, 
if A12 = A21 = 0, then ai = 2 and 02 — 0, thus 
Hc2{T)4:(Tc - T)<j>o/7r'^Di is determined by Di for the 
band with the highest coupling constant An. However, 
if all constants Xmm' are finite, but Di and D2 are very 
different, then Hc2 is determined by the maximum diffu- 
sivity, for example, by D2 ^ Di if the scattering in the 
"weak" band 2 dominates. 

The zero-temperature Hc2{0) can be obtained using 
the asymptotic behavior of U{x) « In^Ajh) for h 00. 
Then all terms proportional to hit cancel out, and Eq. 
( ^ ) reduces to a quadratic equation for In Hc2 , whence 



Hc2iO) = 



<t>oTc ,9-. 



.9 = 



4+ln^?^ + ^ln?i 

w 



1/2 



Xo 

w 



(33) 
(34) 



For Di = D2, Eqs. (g3D-(^ again yield tke result 
Hc2{0) = (l)oTc/2jD of the one-band dirty limitES. How- 
ever for Di ^ D2, Eqs. (|3|)-(34) predict a significant en- 
hancement of i?c2(0) as compared to i?c2(0) = 0.69iJ^2Tc 
for the symmetric case Di = D2. Indeed, in the limit of 
very different diffusivities, we obtain 



Hc2{Q) 



-e 2™ . 



27D2 



D2 < Die- 



(35) 



Hc2iQ) = T^rf e^^, Di « D2e--. (36) 
27£»i 
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FIG. 1: Temperature dependences Hc2{T) calculated from 
Eq. ( ^l| ) for different ratios of D2/D1 and the coupling con- 
stants X^^ ^ 0.81, A™ ~ 0.285, Xa^ ^ 0.119, X^^ « 0.09 
calculated for M gB2 



Unlike Hc2{T) at T « T^, the limiting value of i?c2(0) is 
determined by the minimum difFusivity. It is the feature 
of a two-band superconductor, which causes both the up- 
ward curvature and the enhancement of 11^2 low T in 
Fig. 1. In the limit D2 ^ -Di, the difFusivity D2 does 
not affect Hc2 practically for all temperatures T < Tc, 



except a narrow region T <^ Tc, but for D2 
upper critical field Hc2{T) diverges at T = 0. 



ANISOTROPY OF Hc2 



0, the 



A. General equations 

To address the anisotropy of Hc2, observed in MgB2, 
we consider a uniaxial crystal in a field H inclined by the 
angle 9 with respect to the c-axis. Then the tensors Dap 

in Eq. (^)-(^) have two equal principal values Dm^ in the 
ab plane and a different value D'f^ along the c-axis. If the 
z-axis is taken along H and the a-axis coincides with the 
y-axis, D'^'^f has the following nonzero components 



jjyy 

JJZZ 
JJXZ 



cos^ + sin^ 0, 
cos^ e + sin^ 0, 
^L'"^ = (A^n^-^L°^)sin0cos6 



In the gauge Ay = Hx^ the Uzadel equations (|22 
take the form 



- D^^(9 f 

-D^'^B f - 



2/3)^^(9 f 



2niHx/<l)offn 
+ 2ujfm = 2A„ 



(37) 
(38) 
(39) 
(40) 

and 



(41) 



The solution of this equation has the form fm{T^) ~ 
exp{ikzz - ik^xD^/D^ + ikyy)fm{x), where fm{x) is 
determined by the following equation 



-D: 



-:r:n+Dyji{2iTHxi<t>^Yf„ 

+kl{D^^-Dlf/D^^)!„ 



2ujf„ 
2A 



(42) 



the prime denotes differpatiation over x, and ky was ab- 
sorbed by the shift of xc3. The upper critical field Hc2 
is the maximum eigenvalue of Eq. (p2|). Because of the 
stability condition D^D'^ > D^^, the term quadratic 
in fcz always decreases the eigenvalues of Eq. ( ^^ , so Hc2 
corresponds to k^ — 0. In this case the lowest eigenfunc- 
tion (fo of Eq. (E^) is 



"27rF' 


1/2 


■ nyy' 


1/4 


4>o 




p\xx 





(43) 



An interesting situation occurs if the ratio D^y/Df^ 
is different for different bands, as characteristic of the 
a and tt bands of MgB2- In this case the solution 
fm{x) oc A,„(x) (X (pi^{x) satisfies Eq. (p3|), but does not 
satisfy the self-consistency condition (Q) , because qi ^ 52 • 
This feature of two-gap superconductors can essentially 
complicate the calculation of Hc2 for inclined fields, as 
compared to the case H||c. 

Hc2{0) for an arbitrary field orientation can be ob- 
tained by expanding the solution of the inhomogeneous 
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Eq. ( |43D in the series of orthogonal normahzed eigen- 
functions (p„(gmx) for the different Landau levels n: 



n=0 



w + (2n + 



Am{,u)ipn{qmu)du, (44) 



= (Z?^-i?J^^)i/2^i//0o, (46) 

where Hn{qx) is the Hermitian polynomialsil. Inserting 
Eq. ( ^5| ) into Eq. (^ and summing up over w, gives iJc2 
as the maximum eigenvalue of a matrix M„„' , for which 
the determinant ||M|| — 0. The matrix M„„' is given by 
(see Appendix C): 



(1 - Aii^^W)J„„- - Ai2A2iF^'^ ^ 



(1) ynsVn'sFs 



(2) 



1 — A22-Ff 



(2) 



(47) 



2,.. ^Jl , a„(2. + l)^^ /Ij^ ^^^^ 



ttT 2^r 



= / 'Pn{qix)ips{q2x)dx 



(49) 



where is the Kronecker symbol. For H||c or for 
equal anisotropy ratios D^f^/D'^f' = / D2'^\ we have 
qi = q2, so Vnn' = ^nn' , and all off-diagonal terms of the 
matrix M vanish. In this case the equation for the max- 
imum eigenvalue Hc2 is simply Mqq — 0, which reduces 
to Eq. (^l|) of the previous paragraph. For arbitrary 
field orientation and general relations between the cou- 
pling constants and the anisotropy ratios Dm^ /Dm^ the 
equation ||M|| = for Hc2 can be solved numerically. In 
this case the matrix M„„' should be first diagonalized to 
MnSnn', and then Hc2 can be found as a root of Mq = 0- 
The equation for Hc2 greatly simplifies for the moder- 
ate anisotropy characteristic of MgB2. To quantify the 
degree of anisotropy, we introduce the asymmetry pa- 
rameter C = [(g| - gi)/(g| + ql)?, that is, 



7D: 



(50) 



constant is still a very good approximation. The remain- 
ing summation over s in Eq. ( [47| ) can then be performed 
exactly, using the sum rule due to the orthogonality of 
the eigenfunctions (p„ (see Appendix C): 



E 



(51) 



Because of the 6nn' in Eq. (^l]), all off-diagonal elements 
of the matrix M vanish, thus the equation i\/oo = for 
Hc2 takes the form 



(1-An^^o^'0(l-A22^^o^ 



\ \ p(l)p(2) 



(52) 



In this equation all qm{9) cancel out, and the anisotropy 

manifests itself only via Vim (0) in the functions Fq™'' de- 
fined by Eqs. (^) and (||). After subtraction Eq. (|2|) 
for Ic — hi{2juj]j/T:Tco), from Eq. (|52|), the latter reduces 
to Eq. (|l|) if the following rescaling = {D^Dvyy/^ 
is made. Therefore, all formulas of the previous section 
can be generalized to the anisotropic case of inclined field 
by replacing Di and D2 with the angular-dependent dif- 
fusivities Di{9) and D2{9) for both hands: 



Dm{e) = [D 



= fn(-)2cos2( 



+ D(^^Dl^lhm'9]'/' (53) 



For example, Eqs. ( plD and (53) give the following angu- 
lar dependence of Hc2{9) near T^- 



Hc2{0) 



Tr^[aiDi{9) + a2D2{9)]' 



(54) 



For Di — D2, Eq. (54) reduces to the well-known result 
of the anisotropic GL theory, H^2 = Ho{Tc - TU'k'^D{9) 
for a one-gap superconductor in the dirty limito. 



B. Temperature dependence of the anisotropy of 

Eqs. ( |3l| ) and ( ^3| ) determine both the temperature 
and angular dependences of Hc2- Here the anisotropy of 
Hc2{9) essentially depends on both T and the in-plane 
diffusivity ratio d'^'^ /d[°'\ For instance, let us consider 



a simpler case of isotropic diffusivity Dj"'' = -^'2'^'' for the 
band 2 and anisotropic diffusivity with d[°^^ ^ d'^^ for 
of 7fc2 of MgB2 near indicates that I?^"^ « D^^' and ^-^^ ^^nd 1, which qualitatively models the 3D tt band 

and the 2D a band of MgB2. Then for d''^^ » D^^'^\ 
the angular dependence of Hc2 is most pronounced at 
Tc were it is determined by the large anisotropic Di{9) 
(see Eq. (^4|), while at lower T the angular dependence 
of Hc2 decreases, as it is mostly affected by the nearly 
isotropic small D2 (see Eqs. (^-(|^). However for 



As shown in the next section, the observed anisotropy 



DJ"^ « AD[''\ C « 1/9. In this case the matrix elements 
Vo,2s OC rapidly decreases with s (see Appendix C). By 



contrast, the factor i^i^-'/(l — A22fs^'')) which varies much 
slower with s, can be replaced by its value at s = and 
taken out of the sum. Indeed, near Tc where -~^ 0, the 
function Fg ~^ \ti(2'^ujj:)/ttT) is independent of s over a 
very wide range of s, so the above procedure becomes ex- 
act. Even for T = 0, the function F oc \u[{2s + l)n^/Tco\ 
varies much weaker as compared to the exponential de- 



(2)n 



cay of Vq^, so the replacement of i^i^V(l ~ ^22Fs^') by a 



^(2)^ 



dJ"' < D^^\ the situation reverses: the upper critical 
field is nearly isotropic at Tc, becoming more anisotropic 
at lower T . These two different types of angular depen- 
dencies of Hc2{9) are shown in Fig. 2. 
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FIG. 2: Angular dependences of Hc2{9)Jqr T" « and T = 
calculated from Eq. jsil ) and Eq. (|33|, respectively with 
D*'^^ = Di(6)) = (D("''cos2 6( + D^''£'*"'sin2e)i''2 for 

cases (a) and (b) correspond to different 
ratios / D^^\ and A„„/ are the same as for Fig. 1. 



The above arguments indicate that the temperature 
dependence of the anisotropy parameter H^}2/H^ can 
be rather sensitive to the ratio D1/D2, which in turn 
is determined by the type of intraband scatterers. In- 
deed, if the intraband scattering is dominated by im- 
purities (Hke C or N) which mostly reduce Di in the 
main cr-band, {Di <^ Di), then, according to Eqs. 

( ^ ) and (|33|), the ratio ^^c2/^c2 minimum at Tc 
and increases as T decreases. This behavior has in- 
deed been .nbsfirved .. an . jnairy_ mostly clean, MgB2 
samplesyi^yyHBEiyym. By contrast, if impu- 
rities (like Al, Mg vacancies) cause the strongest intra- 
band scattering in the "weak" 7r-band, {D2 <C i^i), then 
the anisotropy parameter hI^^/H^ would decrease as T 

decreases, The temperature dependence of Hc2/-ffc2 foi' 
these limiting cases is shown in Fig. 3. 




T/T^ 

FIG. 3: Temperature dependences of the anisotropy ratio of 
Hl2/Ht2 calculated from Eq. @ with D2 = Da"' and Di = 
for H^2, and D2 = -D^"' = D[^^ and Di = /d^^ 
with _D*"' = 4:D[^^ for hI^. The cases (a), (b) and (c) corre- 
spond to different ratios Dj"'/^!"') and \mm' are the same 
as for Fig. 1. 
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VI. DISCUSSION 

The results of this work show that the magnetic behav- 
ior of two-gap superconductors can essentiaUy depend on 
the intraband diffusivity ratio r] = /d''^\ This fact 
can have important consequences for MgB2 for which 
Hc2 can be increased to a much greater extent than in 
one-band superconductors not just by the usual increase 
of pn , but also by tuning the ratio of intraband scattering 
rates via selective atomic substitutions on both Mg and 
B sites. Actually, this possibility is naturally provided by 
the electron structure of MgB2 for which substitutions 
on Mg site mostly disturb the out-of-plane Pz orbitals of 
B thus increasing the tt scattering. Likewise, atomic sub- 
stitutions on the B sites mostly increase scattering in the 
2D a band. As shown above, the upward curvature of 

Hc2{T) for Dj""* ^ -Di"'' can make the zero-temperature 
-ffc2(0) significantly higher than what the one-gap esti- 
mate (P suggests. For example, if Z?7r = 0.05-Dcr as in 
Fig. la, the field 7?c2(0) ~ 2.3TcH'^2 could approach 
the paramagnetic limit of MgB2, — TOT for the values 
~ IT/K and 30K observed on dirty c-axis 

oriented MgB2 filmaH3. In fact, recent high-field mea- 
surements on this film gave fa-33.5T and H^^ ~ 48r 
which corresponds ~ 0.12DJ2J. Thus, dirty MgB2 
can exhibit upper criticaLfields significantly higher than 
Hc2{0) ~ 30T C|£-|/V&2S'nL^, which can be very important 
for applicationaid. 

The intraband diffusivities Di and D2 can be either 
calculated from first principles using Eqs. (A22), or ex- 
tracted from measurements of Hc2 and p„. For negligible 
interband scattering, Di and D2 can be obtained using 
Eq. (^) for and the relation between intraband dif- 
fusivities and the residual conductivity a = 1/pn- 



For negligible interband coupling 02 0, Eq. (56) re- 
duces to the second Eq. (l|), of the one-gap theory, 
whereas the stability condition (^) requires that B'^2 
cannot be smaller than AeksNipn/TT. This condition 
gives B'^2 > 2r^T/K for p„ = Imflcm, and Ni = 0.3 
states/cell evH, thus 5^2(0) for Tc = 30X25 cannot be 
smaller than Q.7B'^2Tc = 12.6T even for the symmetric 
case Di = D2 which gives the minimum -ffc2(0) for a 
given H'^2- Foi' asymmetric case, I?^ '^rfi'^ character- 
istic of dirty c-axis oriented Mgiig, filmstj, i?c2(0) can 
considerably exceed this minimurrH. 

Another property for which the diffusivity ratio 

-Dg^V-^i"'' is essential is the temperature dependence of 
the anisotropy of Hc2- As shown in the previous sec- 
tion, the anisotropy parameter H^}^/ increases as T 
decreases for D„ <C D^, , but decreases as T decreases for 
D„ ^ (see Fig. 3). Both cases are in stark con- 
trast with the one-gap theory in which the anisotropy 
ratio for both the upper critical field ff^jj/ffj^ and the 
lower critical field hI^/H^ is determined by the ratio of 
the effective masses rric/mab and is temperature indepen- 
dent. By contrast, the two-gap MgB2 exhibits a strong 
temperature dependence of iJ^^jj/T?, 
ferent from the temperature dependence of H]}-^/. 
was shown both experimentally and theoreticallytHlH 
For the dirtjt, c-axis oriented films, |the observed H^^2l 
is about 2E3, and 13 tt — O.lDo-EJ. In this case the 
anisotropy of Hc2 near Tc is entirely determined by the 
anisotropy of Dg-, thus from Eq. (Q) it follows that 



which is rather dif- 
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Tap 



EN D"'' 



(55) 



For H||c, solving two linear equations ( |32| ) and (55) ex- 
presses Di and D2 in terms of observed i/^2 and p„: 



Di = 

D2 = 



1 



/ 8A:b0oA^2 02 



e'^Pn 



aiN2 - a2iVi V ■^^fiH'c2 

1 f ai _ SkB(l)oNi 

aiN2 - a2Ni 



e^Pn 



n^hH'^2 



(56) 
(57) 



where the normal units are restored. . .. . 

Using the results of ab-initio calculationsoEZl for 
M5B2, \aa - 0.81, A^^ « 0.285, « 0.119, « 
0.09, N„ « 0.3 states/cell eV, and « 0.41 states/ceh 
eV, we obtain A_ = X^cr ~ Atttt — 0.525, Aq = (A^ + 
4A^„A^^)i/2 = 0.564, thus ai = 1 + A_/Ao = 1.93, 
a2 = 1 — A_/Ao = 0.07, and aiN2 » a2iVi. Since both 
Di and D2 must be positive, the slope 3^2 for a given p„ 
in confined within the following limits 
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grateful to D. C. Larbalestier for illuminating discussions. 



APPENDIX A: DERIVATION OF THE UZADEL 
EQUATIONS 

We derive the Uzadel equations for an anisotropic two- 
gap superconductor from the general Eilcnberger equa- 
tions for the quasiclassical Green functions /(v, r,cj), 
/+(v,r,u;) and g{-v,r,uj)B 



n, 



I 



(2w-f vkn)/(k) 

|ukqn5(k)/(q) 



2A5(k) = 
/(k)5(q)], (Al) 



(2a;~vkn*)/+(k) - 2A*g(k) = 



I«kqn5(k)/+(q) - /+(k)5(q)]. (A2) 



8ekB 8ekB 

-Nipn < B'c2 < -N2Pn 



irai 



■na2 



(58) 



Here dAk/wk means integration over the FS with the lo- 
cal density of states 1/wk, the integral terms account for 



9 



nonmagnetic impurity scattering, Uk,q is the scattering 
amplitude, is the density of impurities, is the nor- 
mal group velocity on an anisotropic FS, the wave vectors 
k and q lie on t he F S, an d th e asterisk means complex 
conjugate. Eqs. (Al) and (A_2) are supplemented by the 
normalization condition 



the equation for the superconducting gap 

A(k,r) = 27rT^ / ^y(k, q)/(q, r, c.), 
and the expression for the current density 

J = -AnTelm^ / — 3-Vqg(q, r, w), 



(A3) 



(A4) 



(A5) 



LU>0 



where y(k, q) is the pairing potential. 

Now we define the Green functions /i and /2 on two 
separate sheets of the FS, and write the c ollis ion integrals 
in the right hand side of Eqs. ( Al) and ( A2) in the form 



Stl — Tli 



dA. 



qi 



?^kiqingi(ki)/i(qi) - /i(ki)gi(qi)] 



^|^.k,q.n5i(ki)/2(q2) - /i(ki)g2(q2)] (A6) 

/d A 
^i^|7.k.q.P[g2(k2)/2(q2) - /2(k2).g2(q2)] 
^'q2 

+n, / ^|^,k,q,n52(k2)/i(qi)-/2(k2)gi(qi)] (A7) 



The first integral term in Eqs. ( A6) and ( A7) describe the 
intraband scattering for which the wave vectors k and q 
lie on the same sheet (1 or 2) of the FS, while the second 
integral term describes the interband scattering for which 
k and q lie on the different FS sheets. 

To derive the equations for /^(r, w), /+(k, w), 
f/m(k, r,a;), and Am(r) in the dirty limit, we follow the 
procedure developed for one-band superconductors, ex- 
panding /(k, r,a;) and (7(k, r,w), in suherical harmonics, 
and keeping only the first dipole tcrmcil 

f,n{k, r, uj) = /„i(r, uj) + (v,„(5f„,.(r, cj)), (A8) 
5m(k, r, Lu) = gm{r, uj) + (v„Jg„(r, uj)), (A9) 

where the band index m runs from 1 to 2, and the 
last terms in the right hand side describe first order 
anisotropic corrections to /(k, r, ui) and ^(k, r, w). These 
corrections are small if QttTc ^ v'^^/D„i, where Dm are 
intraband diffusivities (see below). The vector correc- 
tions 6grn and Sfm are related through the normalization 
condition (A3): 



(AlO) 



Now we insert Eqs. ( |A^ ) and ( |A9| ) into Eqs. (|A1|)- 
( |A2| ) and (|A(j )-( |A7|) and integrate over the respective FS 
sheets: / dAk/wm(k), using J Vm(k)dAk = 0. Since all 
/,„ and (5f„i are independent of v, this yields 

2ujh - 2Ai.9i + rfnjf^ = (.91/2 - 32/1)712, (All) 
2u;/2 - 2A2.g2 + rfxijf^ = {g^h - 31/2)721, (A12) 

dAk 



«™(k)^;^(k) 



(A13) 



where Nm = /rf^k/i'm(k) are partial densiti^ of states, 
and 7mm' are the interband scattering rateal3 



dAq dAk 

Nm J Vm{q) Wm'(k) 



(A14) 



Then we substitute Eqs. ( [A^ ) and (^) into Eqs. (|Al|)- 
dA^ ) and (|A6|)-([A7|), multiply them by Vm(k) and inte- 
grate over the respective FS sheets, using the fact that 
/ Vi(k)v2(k)(iylk = 0. In the dirty limit, this yields the 
following relations 

r,f (/™<55f„ - gm5fi) = r«f n^/m, (A15) 

^'^JUl^aL - 9m5f+P) ^ -C^n^/^, (A16) 

Here the tensor TfJ^ is defined by 

TZ? = ^ f l"q.kPC(k)K(k)-^;/^,(q)] ^^^^^^^^ 

(A17) 



N„ 



Wm(k)Wm(q) 



Eqs. (A15) and (A16) can be rewritten as 



fmSgZ - gmSfZ = F^^'^pfm. (A18) 
- 9mSf+" = -F^fuy+, (A19) 



where f,^^ ^ (r^7)-V7/3. Multiplying Eq. (AI8) W 



and Eq. ( |A19| ) by /m, adding, and then using Eq. ( [AlC| ), 
we obtain 



if^^pfm - /™n;/+) 



(A20) 



Using Eq. ( |A15D , ( LA16| ), ( |A20| ), and the normalization 
condition g'^ + |/m| = 1, obtain 



Sfm = ~Fr/igmTlpfm - fmyt^gm). 



(A21) 



Inserting Eq. (A21) and its conjugate into Eqs. (All) 
and ( A12[ ), we arrive at Eqs. ^ and (|) in which the 
intraband diffusivity tensors are given by 



m m \ m J m " 



(A22) 



Averaging the gap equation (A4) over the FS, yields Eq. 
(0) in which the coupling constants Amm' are given by 



A,] 



If dAkdAq 

Nm J t;m(k)Wm'(q) 

whence iViAi2 = Af2A2i. 



^(k,q), 



(A23) 
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APPENDIX B: EQUATIONS FOR Hc2 AND 
WITH THE ACCOUNT OF INTERBAND 
SCATTERING 

For H\\c, the linearized Eqs. (|^) and (^) take the form 



^/l-^n2/i=Ai + (/2-/i)7i2, 
w/2 - = A2 + (/l - /2)721, 



(Bl) 
(B2) 



In the gauge Ay = Hx, the solution of Eqs. (Bl) and 
(H) is 



= /™e-«'"'/2^ A„(x) = A„e-«'-'/2. (B3) 

Here = 2TrH/(j)o, and the amplitudes /,„ can be ex- 
pressed via Ai and A2 from Eqs. (Bl)-(| 



/i = 

/2 



ff+Ai - gi2A2 ^ g_Ai +ffi2A2 
g_A2-g2iAi g+A2-l-52iAi 



a; + ri+ 



(B4) 
(B5) 



where 



20+ = tj+ +7+ + f^o, 2f}_ = tj+ +7+ - f^o, (B6) 



(B7) 



17o = (t^^ +7+ + 27-^-)^^^ 
7± = 712 ± 721, = (-Di ± D2)7tH/4)o, (B8) 

0^- + 7~ 



2g± = 1 ± 



712 721 /T^r,^ 

512 = 7^, 521 = (B9) 

"0 "0 



Inserting Eqs. (B4) and (B5) into Eq. we first sum up 
over u, and express the result in terms of the functions 



F± = 27rr V = In 



^ \2 2ttT J ^ \2 
Then the self-consistency equations (||) reduce to 

GnAi + Gi2A2 =0, 
G21A1 + G21A2 = 0, 

where the matrix Gmm' is given by 



(BIO) 



(Bll) 
(B12) 



Gn = 1 - Xn{9+F+ + g^F^) + \i292i{F+ -~ F_),(B13) 
G12 = Aiigi2(F+ - F_) - Ai2(5-F+ + g+F_), (B14) 

G22 = 1 - A22(g+F+ + g^F^) + A2i5i2(F+ - f^_),(B15) 
G21 = A22.92i(-F+ - F^) - \2i{g+F++g^F^). (B16) 



Using Eqs. (B13)-(B16), the final equation G11G22 — 
G12G21 ~ for Hc2 can be presented in the form 

1 - \ii{g+F+ + g^F^) - \22{g+F+ + g-F_) 
+w{g+F+ + g^F^){g+F+ + g^F^) 
+ (Ai2g2i + A2i.gi2)(F+ - F_) = 0, (B17) 



where Hc2 is th e ma ximum root of Eq. (B17). Setting 
F[c2 = in Eq. (B17), gives the equation for Tc with the 
account of interband scattering for which u)± = = 0, 

^+ = 712 + 721, g+ = 2712/(712 + 721), g- = 2 721/( 712 + 

721). Subtracting Eq. (|l|) for T^o from Eq. Ql?! ), the 
equation for Tc can be reduced to 



U 



P 



712+ 721 \ (Ao-wln<)lni 



27rT J p — w\nt 

712 (Ai - 2A21) + 721 (A2 - 2A12) 
2(712+721) 



(B18) 
(B19) 



where t — Tc/Tco, Ai = Aq + A_, and A2 = Aq — A_. 
The effect of intephaiiA-scattering on Tc suppression was 
considered beforeESSEa 



APPENDIX C: MATRIX M 

Inserting Eq. (^ ) into the gap equation (^ yields 

A„i 2ttT ^ ^ \mn' ^ , . ? /J^ I 1 To ' (^1) 



n=0 



where gI™' = (p„(gma;) Am(x)(ix. Now we multiply 
Eq. ( |Cl| ) by ips{qmx), integrate over x using the orthogo- 
nality condition J ipn{qmx)ipn'{qmx)dx — Snn' , and then 
sum up over lo. This gives two sets of linear equations 



(l-AnFW)GW=Ar2E^-^i'^^^^ 

n 

(l-A22Ff))Gf =A2i^K.F«gW, 

n 

where and F^"^-' are defined as follows 

/oo 
(ps{qix)(pn{q2x)dx 
-00 

pirn, = Y: 

Eq. " 



UJ>0 



1 

uj + {2n+ l)f2ri 



(C2) 
(C3) 

(C4) 
(C5) 



f) can be expressed in terms of the ip function as 
in Eq. (M). Inserting G,l^^ from Eq. (|c|) into Eq. (O) 



^(1) 



yields the closed set of linear equations for Cn 

00 

J2 Mnn'CiV = 0, (C6) 
n'=Q 

where the matrix M„„' is given by Eq. (^^. The solv- 
ability condition of Eq. ( |C6| ) gives the equation for Hc2, 
which is the maximum eigenvalue of the matrix M. 

A general formula for Vnn' is rather cumbersome, so to 
illustrate general behavior of Vnn' we considjCr the case 
of n' = for which the integral ( C4 ) equalalS 



VnO = 



2(?i(72(2n)! 



22"(n!)2(q2 + g2) 



1/2/2 2 

ql + ql 



(C7) 
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if n = 2s, and Vno = for n = 2s + 1, where s is any 
positive integer. As follows from Eq. (C7), V2s,o 
exponentially decays with s. The matrix elements Vnn' 
obey the useful sum rule '^'^\VnsVn's — ^rm' , which can 
be proven by inserting Eq. ( |C4| ) under the sum: 

/ dxdx'(pn{qix)ifs{q2x}'fn'{qix')'fs{q2x') (C8) 

s=qJ-co 

To sum up over s, we use the identity@ 



^ips{q2x)(ps{q2x') 



S{x- 



(C9) 



After inserting Eq. ( |C9| ) into Eq. (|C8|) , the latter reduces 
to / (fn{qix)(pn'{qix)dx, giving the sum rule ([5l|). 

Another equation for Hc2 valid for weak interband cou- 
pling A12A21 ^ A11A22 and any anisotropy can be ob- 
tained as follows. Since off-diagonal terms of M in Eq. 
( ^ ) are proportional to the small factor A12A21, the ma- 
trix element Mq can be obtained in the second order per- 
turbation theory: 



Mn = Ml 



00 



s=l 



M- 



2s, 2s 



(CIO) 



Here Mo,2s is proportional to the small parameter 
A12A21CJ so the sum in Eq. ( CIO ) is by the factor 
^ Ai2A2i^ smaller than the sum in the right hand side 
of Eq. (I3) for Moo. Thus, the contribution of the off- 
diagonal terms to Afo can be neglected, and the equation 
for Hc2 is again A/qq = 0, that is. 



1 — AhFq^-* — A12A21FQ 



r(2) 



1 — A22i^s 



(Cll) 



For moderate anisotropy (C ^ 1), Eq- ( Cll ) redu ces to 
Eq. (p2|), as the convergence of the sum in Eq. (Cll) 
is due to the rapidly decreasing factor V^^ cx (^^ . How- 
ever, for very strong anisotropy ^ « 1 at low T when the 

f 2) I 1 

dependence of Fs on s is essential, Eq. ( Cll ) may be 

simpler for numerical solution than the full Eq. (C6). 
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